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anisotropic superconductor with nonmagnetic impurities is calculated within 
a weak-coupling mean field approximation. It is shown that its dependence 
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There now exists a considerable experimental evidence supporting (i-wave superconduc- 
tivity in the cuprates, but the most direct probes of the superconducting state like the 
electromagnetic penetration depth, photoemission and quantum phase interference measure- 
ments neither confirm nor exclude a possible small s-wave admixture in a predominantly 
dx2_y2 superconductor. tl The linear temperature dependence of the penetration depth at 
low temperaturesB observed in YBCO agrees with the theoretical predictions for a t/-wave 
state.H However, the measurements only went down to about IK and an exponential behav- 
ior below this temperature, indicating a small nonzero gap minimum cannot be eliminated. 
Even by taking data at much lower temperatures, the presence of a small s-wave component 
in the order parameter cannot be entirely excluded in the penetration depth experiments. 
Similar constraints limit the angle resolved photoemission spectroscopy (ARPES) method. 
Although ARPES dataau'S are consistent with a d^z^yi scenario in BSCCO as well as in 
YBCO, the experiments cannot decide with an accuracy greater than an instrumental reso- 
lution if the order parameter completely vanishes at the d^^^y^ nodal lines. This leaves the 
possibility of a small (< 2meV) s-wave admixture. Therefore the above experimental meth- 
ods do not rule out the presence of a small isotropic component, but place an upper bound 
on the minimum of the gap function. As analyzed in Ref. 1, the emerging picture from the 
Josephson experiments supports a scenario of a real mixture of s and dx^-y-i states in YBCO, 
but also does not definitely confirm the presence of the s-wave component. The existence 
of even a small s-wave admixture in a (i-wave superconductor may be tested by thermody- 
namic measurements in the presence of nonmagnetic impurities. It is well known that the 
(i-wave state is strongly suppressed by the defects,EH3 and the s-wave state is not affected 
by the nonmagnetic scatterers.Q In the case of (d + s)-wave superconductor a power-law Tc 
suppression^ should be observed above certain impurity-doping level and the thermodynamic 
properties at large impurity concentration should resemble those of the s-wave state. In fact 
the critical temperature of YBCO is decreased below 12K by the electron-irradiation,EJ and 
the Pr-doping or ion-beam damage lead to a long tail Tc suppression^ characteristic for a 
small nonzero value of the gap function integrated over the Fermi surface.^ However, despite 



the electron-irradiation removing the planar oxygens produces the nonmagnetic defects,llj it 
has not been determined, whether the scattering centers created by Pr-doping and ion-beam 
damage in YBCO are purely nonmagnetic. 

In the present paper we suggest that more significant features attributed to the s-wave 
part of the order parameter may be seen in the specific heat measurements. We study a 
nonmagnetic impurity effect on the specific heat jump at a superconducting-normal phase 
transition in anisotropic superconductors and show that the result depends on the Fermi 
surface (FS) averages of the first four powers of the superconducting order parameter. A 
particularly large difference in the specific heat jump between the states with a nonzero and 
a zero value of the order parameter FS average is observed. We suggest that this measure- 
ment may be used as a test for the presence of an s-wave admixture in a d^-^-yi state. We 
take h = kB = I throughout the paper. 

We consider the effect of potential scattering by nonmagnetic, noninteracting impurities 
on the order parameter with its orbital part defined as follows 

A (k) = Ae (k) (1) 

where e (k) is a momentum-dependent function normalized by taking its average value over 
the Fermi surface (e^) = JpgdSkn{\i.) e^ (k) = 1, where JpsdS^ represents the integration 
over the Fermi surface and n (k) is the angle resolved FS density of states, which obeys 
JpgdSkn (h) = 1. This normalization gives A the meaning of the absolute magnitude of 
the order parameter. The function e (k) may belong to a one dimensional (ID) irreducible 
representation of the crystal point group or may be given by a linear combination of the 
basis functions of different ID representations. The impurity effect is studied in the t-matrix 
approximation.li^GJ This approach introduces two parameters describing the scattering pro- 
cess: c = l/(7rA^o^) and F = ni/nNo, where Nq, Vi and Ui are respectively the overall 
density of states at the Fermi level, the impurity (defect) potential and the impurity con- 
centration. We assume an s-wave scattering by the impurities, that is Vi does not have an 



internal momentum- dependence. It is particularly convenient to think of c as a measure of 
the scattering strength, with c — ;> in the unitary limit and c ^ 1 for weak scattering i.e. 
the Born limit. 

The amplitude of the order parameter is determined by the mean-field self-consistent 
equation 

A(k) = -TEE^(k,kO ,. .tflArvOP (2) 

0. k' ^ +6' + |A(k) r 

where T is the temperature, ^k is the quasiparticle energy, u = 7rT(2n + 1) (n is an integer), 
and V (k, k') is the phenomenological pair potential taken as 

V^(k,k') = -K.e(k)e(k') (3) 

We have assumed a particle-hole symmetry of a quasiparticle spectrum. The renormalized 
Matsubara frequency uJ (k) and the renormalized order parameter A (k) are then given by 

u = uj-Eo, A (k) = A (k) + Si (4) 

with the self-energies defined as follows 
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and the go, gi functions determined by the self-consistent equations 

= ^ V A(k) , . 

To proceed further, we restrict the wave vectors of the electron self-energy and pairing 
potential to the Fermi surface and replace J2k by A^o Jfs dSkfi (k) / dC,k- Integrated over ^k 
the gap equation (|]) can be transformed after a standard procedureli^ into 



where the f^j function is defined as follows 



fu;= dSkU (k) 
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We expand Eq. (H) in powers of A^ around A = using the relations (^-®. Keeping up 
to the fourth power terms in A we get the gap equation in the Ginzburg-Landau regime 
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where the coefficients are given by 
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Taking the derivatives with respect to A^ 
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and with a use of the relations given in Eqs. (|l)-(0) we calculate /o and /i coefficients 
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where g = [r/(c^ + 1)] / (27rT) and ip, tp^^'^ {n = 1,2,3) are the polygamma functions.ll^ 
In the unitary hmit c —>■ and g = r/(27rT). Alternatively for weak scattering (c ^ 1) 
we keep only the terms linear in 1/c^ in a Taylor's expansion which leads to the Born 
approximation scattering rate g = nNoUiV^^ / {2ttT) and g/{c^ + 1) = 0. Coefficients /o and 
/i involve three different types of the Fermi surface averages of the order parameter namely, 
(e), (e^), and (e^). These averages enter the free energy and determine the thermodynamic 
properties at the phase transition. In this paper we discuss a specific heat jump at T^, 
AC(Tc) = Cs(Tc) — Cn{Tc), where CsiT^) and Cn{Tc), respectively are the specific heat of 
the superconducting and normal state, Cn{Tc) = (27r^/3)A^o^c- We obtainM from Eq. (|10|) 
that 
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where ft is ^ at T = T^. This rather cumbersome formula, when considered along with Eq. 



(p!6|), reduces significantly for (e) = case 
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where /i = (1 — c^)/(l + c^). For an appropriate choice of (e^) value, AC {Tc) /Cn (Tc) from 
Eq. (|T9|) agrees with the result obtained by Hirschfeld et al.Ej as well as that by Suzumura 
and Schulzll3 in the Born limit. 

It is informative to discuss the limiting cases of Eq. (|I8|) , that is a pure system where 
gc = and a highly impure one with ft ^ oo in which Tc — > suppressed by the impurities. 
Using a series representation of /i functionE3 we get in gc = limit 

fAC{Tc)\ _ 24 1.426 



(20) 



The Qc ^ oo limit is obtained with a use of Eq. ([I6D and asymptotic forms of polygamma 
functions.El There are two cases to distinguish here. First, when the Fermi surface average 
of the order parameter (e) 7^ then 

/AC(TJ\ 24 
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1.426 (21) 



and the second,li3 with (e) = 0, which leads to 
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We note, that a specific heat jump value in ^c — ^ 00 limit for a nonzero value of (e) given 
by Eq. (^) agrees with that of an isotropic s-wave superconductor. This fact has a simple 
intuitive interpretation. A nonzero Fermi surface average of the order parameter leads to 
an asymptotic power-law critical temperature suppression for large impurity concentration^ 
Tc ~ (Tc„) [F/ (c^ + 1)] , therefore T^ is almost constant for large F values. The 

impurity effect, then, in the large impurity concentration range is the same as in the case 
of s-wave superconductivity, where T^ is not changed by the nonmagnetic impurities. In- 
deed, as it has been shown for the representative order parameters,li2lE2l the gap anisotropy 
is smeared out by the isotropic scattering when (e) 7^ and the density of states approaches 
that of an isotropic s-wave superconductor. Alternatively, for (e) = we observe a strong 
impurity-induced suppression of the critical temperatureBQ leading to a zero value at finite 
impurity concentration, which is reflected by a zero specific heat jump limit value in Eq. 



(22). As a nonzero value of (e) can be achieved only when e (k) contains a component 
belonging to an identity representation of the crystal point group, the measurement of the 
specific heat jump at the phase transition in the limit of Tc ^ (and large impurity con- 
centration) may be used as a stringent test for the occurrence of even a small Aig admixture 
to the order parameter. It should be noted that the effect at large impurity concentration 
for (e) 7^ (Eq. (^)) is independent of the amount of the s-wave content in the order 
parameter, however, as we discuss below, it may be hard to detect for a very small s-wave 



component as it would require an experiment at low temperatures. 

We discuss our results in a context of high- Tc superconductivity, considering a d^2_y2 
state ,EI that is the order parameter given by Eq. (|1|) with e (k) = ik'^ — ky) ( (k^ — ky) ) 
As we have mentioned above, our main result that is the value of the specific heat jump at 
Tc —*> is independent of the amplitude of the s-wave component and its origin. However, 
in order to establish a quantitative behavior of AC(Tc) in a whole range of impurity-doping 
we must work with a particular level of s-wave admixture. We do this by assuming that the 
s-wave component is an artifact of an orthorhombic anisotropy of the system and relate the 
amount of the s-wave admixture to the degree of this anisotropy.BEiJ This approach gives a 
semi-microscopic justification for the {d + s)-wave state. The orthorhombicity in the case 
of YBCO means that the a— and b— crystal axes in the Cu02 planes become inequivalent, 
which leads, with a simple approximation of an elliptical Fermi surface, to the following 
form of an energy banda 

^k = c^kl + Cykl - Ef (23) 

where a ratio of the effective masses Cx/cy is a dimensionless parameter describing the 
orthorhombic anisotropy of the Fermi surface and ep is the Fermi energy. It is easy to see 
within this model, that a {d^^-yi + s) state emerges from a dx2_y2 in a natural way due to 
the orthorhombic distortion of the crystal lattice. A straightforward calculation based on 
a transformation from an elliptical FS to a circular one shows that the normalized drc2_y2 
order parameter defined on the FS given by Eq. ( [23| ) can be represented on a circular FS as 
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where y? is the polar angle. In order to clarify the terminology, we will refer to the circular 
Fermi surface when classifying the superconducting states. Therefore, as a d,j.2_y2 we define a 
state with e (k) proportional to cos2ip and the states with a nonzero s-wave contribution are 
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called {dx2_y2 +s). We note, that the order parameter from Eq. (p^) is d^^-y-i when Cx/cy = 1 
only, that is for a tetragonal symmetry, otherwise it contains a nonzero s-wave component 
proportional to (1 — Cx/cy). In Tab. I we present as the functions of the orthorhombic 
anisotropy parameter Cx/cy the Fermi surface averages which enter the Ginzburg-Landau 
coefficients /o and /i given in Eqs. (|I3|) and (|TBp. We emphasize, that the assumption of the 
orthorhombic asymmetry as the mechanism producing the s-wave admixture in the order 
parameter does not affect the results since only the amplitude of this component matters 
in the calculation. Thus one can obtain the same results in a more phenomenological way 
assuming the presence of the s-wave phase and taking its level as given by (e) in Tab. I for 
the Cx/cy values considered in this paper .Ej 

Based on the discussion of the specific heat jump for a large impurity concentration in 



Eqs. (p])-(p2D we can discuss this limit for d- and {d + s)-wave superconductors. For a 
pure dx2-y2 state {cx/cy = 1) we have (e) = 0. Therefore the specific heat jump decreases 
to zero with a critical temperature driven to zero by impurities as in Eq. (|22D. On the 
other hand, for even a slight s-wave component, (e) ^ and the specific heat jump increases 
and reaches a finite nonzero value at Tc ^ given by Eq. (0). Below, we present the 
specific heat jump at the phase transition normalized by the specific heat in a normal 
state as a function of the normalized impurity scattering rate QcTc/Tcq in the Born limit 
(Fig. la), where QcTc/Tf^ = TrNoUiV^ / {27cTcq) and in the unitary limit (Fig. lb) with 
QcTc/TcQ = r/(27rTc(,). Note that A^o = {cxCy)^^^'^S/{27rh'^), where S* is a sample surface 
area, hence T^^ is different for different values of CxCy product. In the Figs. 2a and 2b we 
show the same AC(Tc)/Cn(Tc) data versus the normalized critical temperature Tc/Tc^. The 
considered states contain a small s-wave admixture varying from about 8% to 16%, therefore 
we observe a strong T^ suppression by the nonmagnetic impurities and a fast decrease in the 
specific heat jump as long as a significant d-wave component is present. Once it is almost 
destroyed and the s-wave part, which is insensitive to the nonmagnetic defects, prevails, the 
BCS normalized specific heat jump value of about 1.426 is restored in a sudden increase of 
AC(Tc)/Cn(Tc). The general tendency of the Tc suppression, given by Eq. ( p!o|) at A = 0, 
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changes at that doping level too and the critical temperature asymptotically goes to zero 
(Fig. 3). For the sake of comparison we show in Fig. 4 the specific heat jump AC(Tc) 
normalized by Cn{Tc) as a function of the impurity scattering rate QcTc/Tc^ in the Born and 
unitary limits for the {s + (i^.2_y2) state, where the s-wave component is large (~ 60%) and 
the dx2-y2 part is considered as minor. 

One can notice from the above figures that the unitary and Born scattering limits differ 
for small values of the pair-breaking parameter and fall on the same curve in the range where 
practically the s-wave superconductivity is only left. The pair-breaking parameter QcTc/Tcq, 
however, has a different meaning in either case. 

We have mentioned before that a detection of a small s-wave component would need a 
measurement at low temperatures. For instance, in a superconductor of the critical tem- 
perature in a clean limit Tc^ = 90K the s-wave content of about 7.4% ((e) ^ 0.074) can be 
observed at a temperature of ~ 2.5K, which is the estimated position of AC{Tc)/Cn{Tc) 
minimum in Figs. 2a and 2b. This minimum is a place where a distinct signal from the 
s-wave component appears, therefore its position is of special interest for possible experi- 
ments. We have found the minimum coordinates (^c^c/^co)* (Fig- 5) and {Tc/Tc^■)* (Fig. 
6) as the functions of the order parameter FS average value (e), which multiplied by 100% 
gives the s-wave fraction in per cent in the normalized to unity order parameter (e^) = 1. A 
plot of (QcTc/Tcq)* vs (e) in Fig. 5 may be also of the experimental use, since the scattering 
rate QcTc/Tcq is proportional to the impurity concentration, which can be estimated in the 
measurements. As one can see from Figs. 5 and 6 the measurements at low temperatures 
are required for small s-wave admixtures, however, they are to be performed at the phase 
transition which should be accessible as long as Tc is measurable. Assuming that a pos- 
sible s-wave admixture is of the order of magnitude of the experimental resolution error 
(~ 2.5 meV) in the ARPES measurements^ of the smallest energy gap values, we can esti- 
mate its fraction as a ratio 2.5 meV/3A meV ~ 0.07, where 34 meV is a measured maximum 
|A| value. Therefore from Fig. 6 we get that the abrupt rise in the normalized specific heat 
jump should be observed at Tc ~ 2.5 K in a superconductor of the critical temperature in 
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the absence of impurities Tc^ = 90K. ExperimentsEj investigating the disorder effect on the 
specific heat jump at T^ in YBCO show AC(Tc)/Tc suppression to zero with the increasing 
impurity concentration. However, the magnetic defects, which act as the pair-breakers on 
both the (i-wave and the s-wave states, were probably present in these studies. 

It is noteworthy that the effect of an abrupt rise in the specific heat jump at T^ may 
be observed even in the purely c?-wave superconductors in the presence of a perpendicular 
magnetic field (H || c-axis). The s-wave component in this case may be induced by the 
vortices.c3 

We have derived the specific heat jump from a mean field weak-coupling theory, neglect- 
ing the fluctuations and the strong-coupling effects. As the observation of thermodynamic 
fluctuations in the specific heat of crystals of YBCO has been reported,E3 we expect our 
BCS result to be modified by the deviations from the mean field approximation. We hope, 
however, that the feature of a sharp upturn in the specific heat will be still present. The 
strong-coupling corrections will rescale the scattering ratesH and may change the magnitude 
of the specific heat jump.E3'EIl 

In conclusion, we have calculated the electronic specific heat difference between super- 
conducting and normal state at the phase transition as a function of nonmagnetic impurity 
scattering rate in the general case of an anisotropic superconductor. We have found that the 
result depends on the symmetry of the order parameter, given by a function e(k), and that 
of the Fermi surface through the following FS averages: (e), (e'^) and (e^). A remarkably 
different dependence of the specific heat jump on the impurity concentration for the systems 
with (e) = and (e) 7^ is observed. We suggest that this effect may be used as a test for 
the s-wave component in the order parameter of the cuprates. 

This work was supported by the Natural Sciences and Engineering Research Council of 
Canada. 
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TABLES 



TABLE L The elliptical Fermi surface averages of the powers of the normalized order parameter 
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FIGURE CAPTIONS 

Fig. 1. Jump in specific heat at Tc normalized by the normal state specific heat at T^ as 
a function of the normalized impurity scattering rate for Cx/Cy = 1 i.e. (e) = (solid), 
c^/cy = 0.9 i.e. (e) ~ 0.0742 (short-dashed), c^/cy = 0.85 i.e. (e) ^ 0.1139 (dot-dashed), 
Cx/cy = 0.8 i.e. (e) ~ 0.1552 (long-dashed): (a) Born limit, (b) unitary limit. 

Fig. 2. Jump in specific heat at Tc normalized by the normal state specific heat at T^ as 
a function of the normalized critical temperature Tc/T^f^ for (e) = (solid), (e) ^ 0.0742 
(short-dashed), (e) ~ 0.1139 (dot-dashed), (e) ~ 0.1552 (long-dashed): (a) Born limit, (b) 
unitary limit. The insets show AC (Tc) /Cn (Tc) in the range of small Tc. 

Fig. 3. Normalized critical temperature Tc/Tc^ as a function of the normalized impurity 
scattering rate for (e) = (solid), (e) ^ 0.0742 (short-dashed), (e) ~ 0.1139 (dot-dashed), 
(e) ^ 0.1552 (long-dashed). 

Fig. 4. Jump in specific heat at Tc normalized by the normal state specific heat at Tc as a 
function of impurity scattering rate for (e) ~ 0.6058 {cx/cy = 0.3) in the Born (dashed) and 
unitary (solid) limits. 

Fig. 5. Position of the minimum in the normalized specific heat jump AC {Tc) /Cn (Tc) on 
QcTc/Tcf) axis (Figs, la, b) as a function of the s-wave component content (e) in the Born 
(dashed) and unitary (solid) limits. 

Fig. 6. Position of the minimum in the normalized specific heat jump AC (Tc) /Cn (Tc) on 
Tc/TcQ axis (Figs. 2a, b) as a function of the s-wave component content (e) in the Born 
(dashed) and unitary (solid) limits. 
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